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Abstract. We present a method of determining a Lax representation for 
similarity reductions of autonomous and non-autonomous partial difference 
equations. This method may be used to obtain Lax representations that are 
general enough to provide the Lax integrability for entire hierarchies of reduc- 
tions. A main result is, as an example of this framework, how we may obtain 
the ij-Painleve equation whose group of Backlund transformations is an affine 
Weyl group of type Eg ' as a similarity reduction of the discrete Schwarzian 
Korteweg-de Vries equation. 



Integrable partial difference equations are discrete time and discrete space ana- 
logues of integrable partial differential equations [H [2], QH [4]- Integrable partial 
difference equations admit classical integrable partial differential equations as con- 
tinuum limits [5j [6l \7\ . Integrable ordinary difference equations are discrete ana- 
logues of integrable ordinary differential equations. Integrable ordinary difference 
equations admit integrable ordinary differential equations as continuum limits [8 . 
Integrable ordinary and partial difference equations possess discrete analogues of 
many of the properties associated to the integrability of their continuous counter- 
parts [gig nsun]. 

We consider partial difference equations whose evolution on a lattice of points, 
wi iTn , is determined by the equation 

(0.1) Q{wi t m,Wl+l,m,Wl,m+l, Wl+l,m+l\ 01,0) = 0, 

where a and j3 are parameters associated with the horizontal and vertical edges 
respectively. The equation is imposed on each square on the space of independent 
variables, (I, m) eZ 2 [3l |4l 1121 1 13] . From a suitable staircase of initial conditions 
|13) . one may determine wi iTn for all (I, m) G Z 2 . Imposing the similarity constraint, 
that 

(0.2) wi +Sum+S2 = wi, m , 

defines a (periodic) similarity reduction |13L 1141 115] . We will assume for simplicity 
that s% and S2 are both positive. In an analogous way to how similarity reductions 
of partial differential equations yield ordinary differential equations [16] , similarity 
reductions given by (|0.2I) yield ordinary difference equations |12l 1131 115] . 
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Given a partial differential equation with some similarity reduction, there is 
a procedure that allows one to obtain a Lax representation of the resulting ordi- 
nary differential equation from the Lax representation of the partial differential 
equation. This holds for autonomous and non-autonomous reductions [17 ]. The 
discrete analogue of this procedure is fairly straightforward for autonomous reduc- 
tions |10[, 111), 118] , however, the process of determining the Lax representation for 
non-autonomous reductions has been somewhat ad hoc |19[, 120] . 

Given a reduction, another task is to determine whether the reduction is a 
known system of difference equations. For autonomous reductions, one may be 
able to find a certain parameterisation which identifies the system as a known QRT 
mapping |21L I22j . which may be classified in terms of elliptic surfaces [23j . For 
nonautonomous reductions, one may be able find a parameterisation of the equation 
that identifies the system as one of the Painleve equations, which are classified by 
the group of symmetries of their surface of initial conditions [24] . 

The aim of this note is to demonstrate a method, which we outline in §1, by 
which we may directly obtain a Lax representation of both autonomous and non- 
autonomous reductions from a Lax representation of partial difference equations in 
an algorithmic manner. The method gives Lax representations in a manner that 
is general and concise enough to directly provide the Lax integrability of entire 
hierarchies of reductions. As an application of this method, we present a reduction 
of the non-autonomous discrete Schwarzian Korteweg-de Vries equation (which is a 
non-autonomous version of Qf =0 in the classification of Adler, Bobenko and Suris 
PUS]) to the g-Painleve equation associated with a surface with A^ip symmetry (or 
q-P(AP)): 

(0.3a) (y'z - 1) [y>z> - 1) ~ 1} ^ " 1} (W " 1} M 1} 



(0.3b) (yz - 1) (y'z - 1) = 



{b iq Hy> - 1) (b 2 qHyi - 1) 
9\ (z — <xi) (z — a 2 ) (z — 03) (z — 04) 
(b\b 2 tz + 6>i) (aia 2 a 3 a4 + OxqHz) 



where t' — q 4 t, the a.;, bi and 9\ are fixed parameters and q is some complex number 
whose modulus is not 1. This is the g-Painleve equation whose group of Backhand 
transformations is an affine Weyl group of type E^p [24] . 

Our method stands in contrast to two methods of performing reductions of par- 
tial difference equations in the literature, namely the method of Hydon et al. [25] 
which is based on the existence of certain Lie point symmetries, and the method of 
Grammaticos and Ramani, who perform autonomous reductions, then deautono- 
mize the equation via singularity confinement [26] . While the first method seems to 
rely on a similar approach to ours, neither method gives rise to the associated linear 
problem for the reduced equation. The approach most similar to our method has 
been discussed by Hay et al. [20 , in which the form of the monodromy matrix for 
autonomous reductions, and its properties, are used as an ansatz for an associated 
linear problem of the non-autonomous reductions of the lattice modified Korteweg- 
de Vries equation. A further extension to this work successfully determined the 
associated linear problem for a hierarchy of systems [19] . 

To demonstrate our method we first provide some simple examples in §2. We 
present an autonomous reduction of the discrete potential Korteweg-de Vries equa- 
tion (dKdV) [7], then present the non-autonomous generalization of this example. 
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= 1. 



mod <?, 



In §3 we first present the g-Painleve equation associated with the surface (oth- 
erwise known as q-Pvi |27| ) as a reduction of (|3.1|) before going to the higher case 
where we present the above-mentioned reduction of (|3.1[) to (10.31) , which we believe 
to be the first known reduction to this equation. 

1. The method 

We start by imposing (|0.2[) as a constraint on our initial conditions, then the 
periodicity gives us that there are si + S2 independent initial conditions to define. 
We solve this periodicity constraint by a specific labelling following |13L 115] ; let 
si = ag and S2 = bg where g = gcd(si,S2), then the direction of the generating 
shift, (c, d), associated with the increment n — > n + 1 is chosen so that 

a 

We specify an n € Z and a p G Z ff by letting 

n = det (; ) p = (Jet [ ' 

W to y 

where the labelling of variables is specified by 

(1.1) Wi tm l-> 

In the case in which <? = 1 the superscript will be omitted. The reduction in the 
autonomous case is a system of g equations given by 

(1.2) Q«,<+t<; C a,<; C a-t«^)=0, p = 0,1,..., 5 -l, 

where a and /3 are constants. In the nonautonomous setting, we have 

(1.3) QK, <+aX+a-t a i. Z 3 ") = °> P = 0, l,...,g-l, 
where a; and /3 m will be, a posteriori, constrained functions of I and to. We will now 
outline how to obtain Lax representations for the autonomous and nonautonomous 
reductions respectively. 

1.1. Autonomous reductions. It is known that multilinear partial differ- 
ence equations that are consistent around a cube are, in a sense, their own Lax pair 
[3l [J- For a generic multilinear equation, (jO.ll) . that is consistent around a cube, 
a Lax pair may be written as 



(1.4a) 
(1.4b) 
where 

(1.5a) Lr m = A; jTO 



(1.5b) M Lm = fii. 



<Pl+l,m — Li tTn (pi tm , 
4>l,m+l = Ml m <j>l m , 



dQ(x, u, v, 0; a, 7) 
dv 

d 2 Q (x 7 u,v,y;a,7) 
dvdy 

( dQ(s,M,u,0;/3,7) 
du 

d 2 Q(x,u 7 v,y;(3,-f) 



V 



dudy 



— Q(x,u, 0,0; a, 7) 
dQ(x,u,0,y;a, 7 ) 

-Q(x,0,v,0;/3 )7 ) 
dy 



X = Wi !rri 
U = Wi + i 



Wl y m+1 
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where 7 is a spectral parameter. The compatibility condition is 
(1-6) Mi + i^ m Li. m — L^ m+ iMi :7n , 

forcing the prefactors, Xi iTn and p,i m , to be chosen in a manner that satisfies the 
equation 

det _ det M; + i ;TO 

detL ijm detM iim 

When the prefactors are appropriately chosen, imposing (ll.6[) is equivalent to (jO.ip . 
In practice, it is often computationally convenient to deal with some transformation 
of this Lax pair. 

To obtain a Lax representation for the system of ordinary difference equations, 
(|1.2[) . we define two operators, A n and B n , associated with the shifts (Z,m) — > 
(I + Si,m + S2) and the generating shift, (I, m) — > (l + c,m + d), respectively. These 
operators have the effect 

(1.7a) 4> n = A n 4> n , 

(l-7b) <t> n +i = B n 4> n , 

where one representatiorQ, that is simple to write, is as follows: 

S2 — 1 Si — 1 

A n <-< J { Mi +Sl . m+1 Yl 

j=0 i=0 
d-1 c-1 
B n <-> Y\_ Ml+c : m+j JJ Li +i ^ mi 

where the dependence on n and p is specified by 

Li, m (wi, m ,wi +1 , m n) i-> ^(7) = ^«,wf^;7), 
M /im (t£>i iro ,«7i + i lW ,;7) ^ MP( 7 ) = A/P«,<; c a ; 7 ). 
The compatibility condition, 

(1.8) A n+ iB n — B n A n = 0, 



is equivalent to imposing (11.21) . We call A n the monodromy matrix for the following 
reason: by identifying all points in Z 2 that are multiples of (si,S2) apart, we may 
consider the space in which the new system exists as being cylindrical. We wrap 
around in a manner that connects the points that are identified by the similarity 
reduction. The monodromy matrix, rather than presenting a trivial action as (|1.7a[) 
suggests, expresses the action of wrapping around the cylinder, as in figure [TJ 

The monodromy matrix can be expressed as a function of the si + S2 initial 
conditions, (w", w° +1 , . . . , wf7;,_ 1 ), by following the standard staircase. Geomet- 
rically, the standard staircase is the path between two lines which squeeze a set of 
squares with the same values, i.e., a set of squares shifted by (si, S2) |15j . 

One advantage of the generating shift is that every other shift in n may be 
expressed as some power of the generating shift by construction [28j. Furthermore, 
this generating shift allows us to constrain the non-linear component, where we 
need to use (|1.2[) . to just g places. We have illustrated the standard staircase and 
generating shift in figure [5J 



4n practice, the product follows the path of a standard staircase 28 . 
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Figure 1. A pictorial representation of the way in which the 
monodromy matrix wraps around to similar points for a (12,4)- 
rcduction. 
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Figure 2. A (9, 6)-reduction and the labelling of variables. In 
this example, the shift (p, n) (p + 1, n) corresponds to the shift 
(a, b) — (3,2) and the shift (p, n) — > (p,n + 1) corresponds to the 
shift (c,d) = (1,1) 



In the example defined by figure [21 if we allow our monodromy matrix to follow 
the standard staircase, the monodromy matrix is 

A n <-Hi; + 8^„ + 6Afi + g !m+ 5Li + 7 im+ 5Li + 6 !m +5Afi + g !m+ 4i; + 5^ m+ 4Afi + 5 !m+ 3 

£z+4,m+3^i+3,m+3A/z+3 im +2£z+2,m+2Mi+2,m+l£; + Um+li(,m+lMi !m , 

and the other half of the Lax pair is 



,m • 
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The resulting compatibility condition, (|1.8[) . gives the evolution equations for the 
K+n i = 0, 1,2: 



QK- 


2i 


K 


-4- 


W n+1' 


W n- 


l'i 


a, (3) 


= 0, 


Q{wl- 


2- 


W° 


-4' 


™l+ll 


«£_ 


li 


a,p) 


= 0, 


Q«- 


2l 


K 


-4' 


™l+ll 


W° 


l'i 


a, /3) 


= 0. 



In general, this procedure gives us an s\ + S2 dimensional mapping, 



which, applied to (u>° , . . 



<\> : £ Sl+S2 

' ! W n+a+b- 



— > 



•"S1+S2 



i), s ives K+iX+2v 



,J n+a+b) 



This 



new set of values forms a new standard staircase. As a matter of fact, this new 
standard staircase is the old one translated by the generating shift. 

1.2. Nonautonomous reductions. To deautonomize this theory, we con- 
sider the a and (3 to be functions of I and m. As L\, m and M/ j?n are shifted in 
only m and I respectively in the compatibility condition, (jl.6|) . replacing a and (3 
with ai and /3 m , which are arbitrary functions of / and m respectively, preserves 
the Lax integrability. Hence, our basic non-autonomous lattice equations may be 
considered to be of the form 

(1-9) Q(wi,m,Wl+l,m,Wl,m+l,Wl+l,m+l](Xl,Pm) = 0, 

where the Lax representation is specified by (|1.4j) where 
(1.10a) 

_ dQ(x,u,v,0;ai,j) 

A/ ' 



dv 



-Q(x,u,Q,Q;ai,~f) 



d 2 Q(x, u, v, y; a u 7) dQ(x, u, 0, y\ a h 7) 



(1.10b) 



dvdy 



( dQ(x,u,v,0; /3 m ,j) 
du 



dy 



-Q(x,0,v,0;/3 m ,j) 



X = Wl, m 

u = Wl+i,m 



V" 



d 2 Q(x, u, v, y; (3 m , 7) <9Q(x, 0, v, y\ (3 m , 7) 



dudy 



dy 



x = Wl >m 

where 7 is a spectral parameter and the prefactors, A/ jjn and ^t/, m , are chosen to 
satisfy the compatibility conditions, in an analogous manner to the autonomous 
case. 

If one assumes that the a and (3 are functions of both I and m, i.e., a — a;, m and 
j3 = f3i, nll then demanding that ai, m is independent of m and (3i, m is independent 
of I has also been shown to be a necessary condition for singularity confinement for 
equations in the ABS list |29| . The above constitutes a Lax pair interpretation of 
this constraint. 

Let us now specialise our choice of systems to those that admit representations 
of the additive form 

Q(wi, m ,wi + i im ,wi tm+ i,wi + i^ m+ i]at - p m ) = 0, 

or the multiplicative form 



0, 
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with a possible additional dependence on ai+i — ai and /3 m +i — Pm in the additive 
case, or oti+i/oti and (3 m +i/Pm in the multiplicative case. A list of transformed 
equations appears in tableQ] where subscripts m and a denote those functions, (jO.ll) . 
dependent on a multiplicative or additive combination of a/ and /3 m respectively. 

ABS Q(x,u,v,y;ai,l3 m ) 

Hl a {Wl, m - Wl+l,m+l)(wi+l,m ~ Wl,m+l) + Pm ~ Oil 

xx -i ( Pm+1 \ ( Pm+1 \ . . O-f 

Hl m Wim t02+l,m+l Wl+X,m 7, Wl,m+1 J I 1 ~7pr~ 

V Pm / V Pm J Pi, 

xx ( Pm+1 \ ( Pm+1 \ «f 

n2 m Wi m Wt+im+i Wi + \ m U>2, m +1 - — 

V Pm J \ Pm J Pi 

+ ( 1 — ■q~ J ( W(,m + Wl+l,m + (wi,m+l + U>Z+l,m+l) J + 1 

\ Pm / \ Pm J 

,m^Z,m+l ^l+l, m^Z+l,m+l) 



H3 m 02 I ^.m^+l.m + o 2 ^,m+l^+l,m+l I 
Pm \ / 



Jm+1 



,m^Z,m+l ^+l,m^+l,?n+l ) "1" 



Ql m -^-(ttfj.m ~ 1"i,m+l)(^+l,m ~ ^J+l.m+l) 
Pm 

-{Wl,m - Wl+X,m){wi,m+l ~ Wl+l,m+l) 



O1 <5^0 a l ( Pm+1 \ ( Pm+1 

Q 1 ™ a w l,m+l Wi+l,m 7, Wj+l ,m+l 

Pm V P™ / V Pm 

Pm+i , w <5af / of \ 

(WZ,m - W(+l,m ) (W*,m+1 - Wi+l,m+l j + "oj" I 1 I 

Pm V Pm / 



A, 



/v , a l ( Pm+1 \ ( Pm+1 
^i 2 m — \ w l+l,m fli—Wl+\,m+\ ) \ Wl, m jjj— Wi,m+1 



P" 2 

■ ™2 (Wl,m - Wl+l,m) (W|,m+1 ~ Wl+l,m+l) 

2 2 

'7-h 1 I u 'i,m + U)i+l,m + 02 ^.ro+l + R2 Wi + l,m+l 

Pm \Prn / \ P m P TO 

m ( on \ f af ai 

M ?2 „ +1 

m Pm 



/3m Vp'm / \/5 m A 



Table 1. A list of various lattice equations (taken from 3, 4 ) in 
a suitable form for non-autonomous reductions. 



With each lattice equation written in terms of on — /3 m or on / (3 m , the necessary 
requirement for (10.21) to be consistent is the requirement that 



Oil — Pm — &l+si ~ Pm+s 

in the additive case, and 

Oil _ Oil+ S1 

Pm P III ■ Sj 
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in the multiplicative case. By a separation of variables argument, we define h and 
q by letting 

(1.11a) ai +si -ai = /3 m+S2 - f3 m := habg, 

(1.11b) ^±£i = ^±£2 . = g «6 9) 

Oil Prn 

in the additive and multiplicative cases respectively. Although it is not a technical 
requirement, we will assume that h is not and that q is not a root of unity. We 
solve the additive and multiplicative case by letting 

ai = Mb + a/, f3 m — hma + b m , 

Oil = a iq bl , j3 m = b m q am , 

where a; and b m are sequences that are periodic of order s\ and S2 respectively (not 
related to the constants, a and b). This choice of ai and /3 m ensures the consistency 
of the reduction with as many degrees of freedom as the sum of the orders of the 
difference equations satisfied by ai and /3 m , (ll.llap and (jl.llbp . i.e., si + S2- 

To provide a non-autonomous Lax pair, we need to choose a spectral variable, 
x, in a manner that couples a linearly independent direction with the spectral 
variable, 7. While any linearly independent direction may be considered a valid 
choice, we present a simple choice. Our choice of spectral parameter is specified by 
introducing the variable k = I and x — hbk — 7 in the additive case and x = q hk /7 
in the multiplicative case. In the additive case 

Ll.rn = Ll,m{m ~ l) ^ L n (di + x) , 

M Lm = M Um {{fi m - ai) + (ai - 7)) ^ M n (x + hn + b m ), 
and in the multiplicative case 

Li, m = Li, m (ai/j) = Li, m (aix), 
M Um = Mi tm ((f3 m / ai )(ai/j)) = M Lm (b m xq n ). 
This gives us a non-standard Lax pair, which, in the additive case reads 

Y n (x + abgh) = A n (x)Y n (x), 
Y n+ i(x + cbh) = B n (x)Y n (x), 
and in the multiplicative case reads 

Y n (q ab 3x) = A n (x)Y n (x), 
Y n+1 (q cb x) =B n (x)Y n (x), 

where 

s 2 — 1 si — 1 

(1.12a) A n (x) 4 i Yl Mi +Sum+j Y[ U+i,m, 

j=0 i=0 

(l-12b) B n {x) 4 1 Yi M i+c,m+j 11 L l+hm . 

j=0 i=0 

The compatibility conditions, 

A n+ i(x + cbh)B n (x) = B n (x + abgh)A n (x), 
A n+1 (q cb x)B n (x) = B n (q ab 3x)A n (x), 
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in the additive and multiplicative cases respectively, gives us (|1.3j) . This choice 
of spectral variable has the advantage that the spectral matrix and deformation 
matrix, A n (x) and B n (x), have a simple dependence on the independent variable, 

71. 

2. Some simple examples 

In this section, we present some examples of the theory above. An example that 
has appeared recently is the example of q-Pvi as a reduction of the discrete modified 
Korteweg-de Vries equation [30 , which also gave rise, via ultradiscretization, to the 
first known Lax representation of u-Pyi- 

2.1. Autonomous example. We consider some additive examples, in partic- 
ular, we will consider reductions of the discrete potential Korteweg-de Vries equa- 
tion, 

(2.1) {Wl,m - W; + l,m+l)(w; + l,m ~ W(,m+l) = OL - /3, 

labelled as Hl a in table [TJ which possesses a Lax representation of the form (|1.4j) 
where L^ m and M;. m are specified by 



(2.2a) 
(2.2b) 



Ll^m — 



M, 



1 



a 







7 - Wl im Wl + l : 
-Wl+l,m 



7 - Wl,mWl^ m +l 
1 -Wl,m + 1 

Let us consider a reduction, (|0.2|) . where si = 2 and s% = 1, with a labelling 
indicated in figure El This gives us g = 1, a = 2 and 6=1, hence n — 2m — I, and 
the direction that characterises the generating shift, (c, d), is chosen to be (1,1). 









1 












W4 


w 3 








w 3 


W 2 






w 3 


Wi 


Wi 


W 






Wi 


w 









Figure 3. The labelling of initial conditions with (2,1) periodicity 
and an evolution in the (1, Indirection. 

The product formula for the monodromy matrix, A„, and the matrix that is 
related to the generating shift, B n , are 

A n =Mi + 2, m Ll+l ym Li ym 

w n P - 7 - w n w n+2 \ f w n+ \ a - 7 - w n w n+ i 

1 -W n+2 J \ 1 -W n 

w n+2 a - 7 - w n+ iw n+2 
1 -w n+ \ 

W n +1 /3 - 7 - Wn+lWn+3 \ ( OL - 7 - W n +lW n +2 

1 -Wn+3 J \ 1 -Wn+1 
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The compatibility condition, given by (|1.8[) . reads 

,m+ ,171-^1, m 

= ^1+3^+1^1+2^+1^1+2 ,m-^'/+l,m-^ y 2,m' 

This simplifies to 

l,m — 

M l+2 

w n +z ol - 7 - w n+2 w n+3 \ / «; n+ i /3 - 7 - u>„ + iw„ +3 



u>„ /3 - 7 - w n w n+2 \ ( w n+ i a - 7 - w„w„+i 

1 -W n +2 J \ 1 ~Wn 

which defines the evolution of this autonomous reduction to be given by the equation 

(2.3) (w n - w n+3 )(w n+ i - w n+2 ) =a-/3. 
If we let y n — w n — w n +i, this equation is equivalent to 

. . a — ft 

(2.4) y„_i + y n + y n+1 = , 

Vn 

which is a well known example of a second order difference equation of QRT type. 

2.2. Nonautonomous example. The autonomous equation and Lax repre- 
sentation generalise naturally to the non-autonomous case by replacing a and j3 by 
ai and j3 m respectively. Furthermore, we may satisfy the periodicity constraint, 

- Oil = Pm+l - Pm ■= 2h. 

We solve this constraint by letting 

ai = hi + ai, P m =2hm + b m , 

where ai is periodic of order two and b m is constant, and hence, may be taken to 
be without loss of generality. The evolution equation for this system may be 
represented as an application of the nonautonomous version of (|2.1j) translated by 
the vector (1,1); 

(2.5) (w n - w n+3 )(w n+ i - w n+2 ) = a l+1 - P m +i. 

recalling that n = 2m — I. The increment in I and to by 1 directly corresponds to 
the increment in n by 1. In the simplest case where ai — a\ is constant (rather 
than periodic), letting 

Vn =W n — W n +1, 

results in the evolution equation 

a i+ i - fi m+ i -hn -h + a\ 
y n + y n +i + y n +2 = = 1 

Un+l Vn+l 

or alternatively 

,„ —hn + a\ 

(2.6) y n _i + y n + y n+ i = . 

Vn 

To form the Lax pair for this reduction, we choose a spectral variable to be 
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Since ai and /3 m appear with 7 in L^ rn and M/, m , this gives us 

on — 7 = hi — 7 + a\ = a\ + x, 
fim — 7 = 2hm — 7 = hn + x. 

which means we may interpret our I and m variables in terms of x and n variables. 
The (I, m) — > (1+2, m+1) shift and (/, m) — > (l+l, m+1) gives us (n, x) — > (n, x+2h) 
and (n, x) — > (n + 1, x + /i), hence, we have a linear system of the form 

Y n (x + 2h)=A n (x)Y n (x), 
Y n+ i(x + h) = B n (x)Y n (x), 

where 

A n (x) ^Mi + 2, m Li + i t7n L Lm 

w n hn + x - w n w n+2 \ / Wn+i a 2 + x - w„w Il+ i 



1 "W„+2 J \ 1 -w, 

W rl+2 dl + X - W n+1 W n+2 
1 -Wn+l 

B n (x) <r-<Mi +lim Li im , 

w n+ i hn + x - w n+1 w n+3 \ / w rl+2 ai + x - w n+1 w n+2 
1 -W„ +3 y ^ 1 -Wn+l 

The compatibility condition is 

(2.7) A n+1 (x + h)B n {x) = B n {x + 2h)A n (x), 

which gives (|2.5[) . To express everything in terms of y n = w n — iu n +i, we use the 
matrices obtained from applying a gauge transformation, L\^ m — ► \ m Li im Si im 
and Mi tm -> S^ +1 Mi !m Si, m , where 



S,. 



1 Wl,r. 

1 



in which case A n (x) and B n (x) are given explicitly in terms of products of matrices 
with entries expressible in terms of the y n variables: 

-Vn-Vn+i (y n + Vn+if + hn + x \ / y n a x + h + y 2 n + x 

l -y n - y n +i J \ 1 2/™ 



A. (a;) 



Bn(s) 



y n+ i 01 + y r 2 l+1 + x 
1 2M+i 

-i/n+i - Un+2 (j/n+1 + 2/n+2) 2 + hn + x 



-Vn+l ~ Vn+2 



y n+ i 01 + y^ +1 + a; 

1 2/n+l 



for which (|2.7p now gives (|2.6p as required. We recover Lax matrices for (|2.4p by 
letting 

A„ := lim i n (i), 
S n := lim B n (x), 



12 



C. M. ORMEROD, P. H. VAN DER KAMP, AND G.R.W. QUISPEL 



B n (x) 



whose compatibility, (|1.8p . gives ()2.4j) where a — /3 = 01. This is a way in which 
the autonomous and nonautonomous reductions are related. 

While the simple case above demonstrates the basic mechanisms in the method, 
the periodicity constraint allows us to build in an extra variable. If we allow the 
full generality of a periodic value of ai , by letting 

f hi + ai where I is odd, 
I hi + a 2 where I is even, 

we obtain 

A n (x) <-lM; +2 ,mi!+l,mii,m 

. . _ / w n hn + x - w n w n+2 \ ( w n+ i a 2 +x - w n w n+1 
An[X) ~ \ 1 -w n+2 ) { 1 

W n +2 ai + X + W n+ iW n+2 
1 -Wn+l 

w n+ i hn + x - w n+ iw n+3 \ ( w n+2 ai+x - w n+1 w n+2 

1 -W n +3 J \ 1 -Wn+l 

The compatibility condition needs to take into account that the shift, (n, x) — ¥ 
(n + 1, x + h), also shifts the position on the lattice, hence, we need to couple 
the non-linear component with a swapping of the roles of a\ and a 2l hence, the 
compatibility condition, given by (I2.7|) . results in the evolution 

(w n - w n+3 )(w n+1 - w n+2 ) = hn + a 2 , 

a i a 2 + h, a 2 — > ai — h. 

If we restrict our attention to even powers of this map, the Oj are constant. We 
proceed to specify a change of variables, which is motivated by some historical con- 
text. In the 1980's, Novikov and Veselov formalized the derivation of the Hamil- 
tonian structure of hierarchies of soliton equations from their Lax representations 
|31L 132] . The Darboux coordinates are the poles of the Baker- Akhiezer function 
and the eigenvalues of the spectral matrix at these poles, which for 2x2 matrices 
coincides with the roots in the spectral variable of the off-diagonal elements and 
the diagonal elements evaluated at those roots. The Hamiltonian description of 
the isomonodromic deformations of [33j are expressed in these coordinates. While 
the link between the symplectic structure for discrete Painleve equations and the 
discrete isomonodromic deformations is not as well developed, it is interesting to 
note that the parameterizations of many known discrete Lax pairs are provided in 
terms of these coordinates [341 1351 127) . Perhaps a discrete analogue of the Baker- 
Akhiezer function holds the key to linking the geometric theory and the theory of 
discrete isomonodromic deformations. It may also provide an algorithmic manner 
of describing the symplectic structure for hierarchies of equations. 

Our new coordinates are (y n , z n ) where y n is the root in the spectral variable 
of the (2, l)-element and z n is a variable which parameterizes the diagonal elements 
at this root. In the lattice variables, w n , these are 

Vn = K+2 - W n )(w n+2 - W n +l) - a 2 , 
_ K+2 - W n )(y n + «l) 

y n + nh 
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from which we extricate the second-order system 

Vn + Vn+2 = %n ~ (Ol + 2 ), 

_ {Vn+2 + ai)(y n +2 +0.2) 
ZnZn+2 - (y n+2 + (n + 2)h) ■ 

This system is equivalent to 

(y n + a 1 ) 2 (y n + a 2 ) 2 



(y n + Vn-2 + ai + a 2 )(y n + Vn+2 + a\ + a 2 ) = 



(y„ + nh)' 2 



which is a special version of cZ-Pjy [36j . The second power of the generating shift 
is equivalent to the shift (I, m) — > (/, m + 1), hence, a simplified Lax representation 
for this system is 

Y n (x + 2h) = A n (x)Y n (x), 
Y n+2 (x) = B n (x)Y n (x), 

where B n {x) -H M^ m . We may simplify the spectral matrix, via a gauge transfor- 
mation, to be 



A n {x) 



(y n + ai)(y n + a 2 ) 2 

a:-2/n (y n +nh)z n 



where 



e =(2/n + /i«)(j/n + «i) + + Vn + a>i)a 2 , 
and a deformation matrix, under the same transformation, becomes 

B n (x) = { 1 

Using this method, one is able to provide product formulas for the Lax represen- 
tations of q-Pn, q-Pv °f Hay et al. [20] and the Lax pair for a version of g-P/jj 
of Joshi et al. [37] . Furthermore, this shows and explains that the Lax pairs for 
the reductions of [20J factorize in a nice way. The resulting factorizations provide 
a simple way to compute the compatibility. 



3. Reductions of dSKdV 

We consider similarity reductions of the nonautonomous discrete Schwarzian 
Korteweg-de Vries equation; 



(3.1) m 



1 1 



1 1 



Wl + i. m +l ~ ^/+l,m w l,m ~ W^ m+ i 
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which possesses a Lax representation of the form (|1.10|) where 
(3.2a) L hm 



1 

oil 



(3.2b) 



Mi, m = 



1 

,m ^/,m+lj 



1 



3.1. q-Y(A^). Recently, one of the authors derived a q-analogue of the sixth 
Painleve equation (or q-P(A^')) as a reduction of the non-autonomous modified 
Korteweg-de Vries equation [30] . This work demonstrated the specified method 
where the transformation between the lattice variables, the u>j, and the Painleve 
variables was relatively simple. We provide a similar but more complicated relation 
between the lattice variables for the discrete Schwarzian Korteweg-de Vries equation 
and the Painleve variables of q-P(A^). 

We now provide a new reduction of the lattice equation to a version of q-Tvi 
given by 



(3.3a) 



(3.3b) 



yy 



b 1 b 2 {z-l)(g 2 z-l) 

1 q 2 {hb 2 t - e lZ ) (hb 2 t - e 2Z y 

(6lgV - 1) (b 2 gHy> - 1) 
q 2 (oitf' ~ 1) M ~ 1) ' 



t' 



q\ 



where Q\Q 2 = a\a 2 b\b 2 . This extends our previous result in [30] by the addition 
of an extra parameter, which represents an integral that is used to reduce the 
order of the map, given by Q\ (or 9 2 ). Obtaining a g-analogue of the sixth Painleve 
equation from a discrete analogue of the Schwarzian KdV equation has some historic 
significance as it could be considered a discrete analogue of the reduction of the 
Schwarzian KdV equation to the sixth Painleve equation |38l 139) . 

We consider a reduction of (|3.1[) . given by (|0.2[) where si = s 2 — 2. The 
labelling, given by (jl.ip . is depicted in figure HJ The constraint, (jl.llbl) . becomes 



(3.4) 



Oil+ 2 fim+2 



Oil 

which introduces the parameter q. 























w° 2 


































w\ ' 
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Figure 4. The reduction and the labelling of variables. 
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We satisfy the periodicity constraints of Q3.4p by explicitly setting 
aiq L if I is even, _ J b\ q m if m is even, 



m ' a 2 q l if I is odd, rm \ b 2 q m if m is odd. 

Naturally, the dependent variable, n, is invariant along the direction of the reduc- 
tion. The correspondence between the associated linear problem here and that of 
Sakai [35j is made more natural by specifying an independent variable, t, and our 
spectral variable, x, following previous sections: 

t = q m - l = q n , x =l = £. 

7 7 

For the sake of clarity, the reduced lattice variables may be regarded as functions 
of t or n under the identification 

w l n+k ^w l (q k t), 

so that the shift n — > n + 1 is equivalent to the shift t — > qt. We know that the shift, 
(I, m) — > (I + 2, m + 2), is equivalent to the shift (x, t) — > (q 2 x, t). Since a — 1 and 
6=1, the direction of our generating shift is (I, m) — > (l,m+ 1), which is equivalent 
to the shift (x,t) — > (x, qt). This gives us a Lax pair of the form 

(3.5a) Y(q 2 x, t) = A(x, t)Y(x, t), 

(3.5b) Y(x, qt) = B(x, t)Y{x, t), 

where the matrix A{x, t) is given by the product 

(3.6) A(x,t) 



(3.7) A(x,t) 



1 / 1 <"<+! 

a 2 x 



w 



n+l 



a\x 




1 - W n \ I ' W n ~ W n+1 ' 



1 



■J n+X ~ w n 

The deformation matrix, B(x,t), corresponds to the shift m — > m + 1, hence, is 
given by 

(3.8) B{x,t)^M l>m , 

( 1 K - v 

(3.9) B(x,t) = txbi 1 

\K -w° +1 
The compatibility condition, written as 

A(x,qt)B(x,t) = B(q 2 x,t)A(x,t), 
yeilds the evolution equation 

,„ , n n _ a ^ W n i W i-l ~ W n) + tb l ( W n ~ W Ll) W n j b 2 

(o-tua,) w n+1 — j-^ 1 \ i j-u t n i — \ — ' 1 — ' 

,„ , nU N 1 a 2 K-l - O + *&2< (K - <_i) 

(3.10b) w n+ i= 7-q ^r— 77 7 i o - S ' b 2 ->qh. 

Notice that the generating shift swaps the roles of b\ and 62, coupled with multi- 
plicative factors introduced to compensate for the dependence of (3 m on q m . We 
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e 2 t 



claim that the second iterate of this mapping is q-P (A^). To make the full corre- 
spondence with (|3.3p . let us expand out (|3.7[) to give a matrix of the form 

A(x, t) = A Q + A x x + A 2 x 2 , 

where Ao — I and the eigenvalues of A2 are 6\t and 6^, where 

9 t = tb ib 2 «+i - w n ) « - wj +1 ) 

( w n - Wn+l) { W n - W r\+l) 
«1«2 K ~ W°+l) (wj ~ wj +1 ) 
t «+!-<)(<-<+,) ' 

where 6i = 8i(n) satisfies 

Oi(n + l)Vi{n) = , 

q 

9 2 (n + l)9 2 {n) = qaia 2 bxb 2 , 

hence, 9i(n + 2) = 9i(n). That is to say that 9\ and 9 2 are 2-integrals of the 
generating shift |40j . with the additional constraint 

6i6 2 = aia 2 b\b 2 . 

Secondly, from the product form, we have that 

det A(x, t) = {a\x — 1) (a 2 x — 1) (bitx — 1) (b 2 tx — 1) . 

In accordance with the motivation given in the previous section, we choose to 
parameterize the spectral matrix, A(x,t), in terms of the x-root of (2, l)-element 
of A(x,t), and the diagonal entries at that root. This provides us with a spectral 
matrix in the curious form 

'tx (x - y n + ( n ) 6>i + 1 txS n uj n 9 2 
(3.11) A{x,t)=\ tx{x~y n )9 x 

1 - tX (X - y n + Tin) 9 2 + 1 



where 



{Vnai - l)(y«a2 - 1) , 1U , , 
, Vn = (bity n - l)(b 2 ty n - l)z ri 



Fixing the determinant requires that 
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Under this identification, the y n and z n are specified by the lattice variables 
(3.12a) 

Vn = - (Wn - W n+1 )((w n+1 - W n )(a 2 (w° +1 - W°) + ht(w° - W n+1 )) 

+ a lK° - W «+l)K° - W n+l)) + - Wn+l)(wl +1 - V)D(wl - W n+1 ) 

~*~ M(™°+1 - w\) (02K 1+ 1 - Wn+l)( w n+l ~ W n) + ~ 



( w n ~ w n+l)) ~ a ia2 (wl[ - W n+1 )(w n - wl +i y 



(3.12b) 



- tb 2 (w° n - w n+1 ){w n+1 - wl) - ht(w n+1 - wi +1 )(w° +1 - w l n )) 

{w° - w° +1 )K° l+1 - w n )((w n - wl +1 )(axw° - a 2 w° n+1 
+ (a 2 - ai)wi +1 ) + b 2 t(w° +1 - - w n+1 )) 



We may now make the correspondence with (|3. 31) via the identification that y = y n , 
y' = y n +2, z = z n and z' — z n+2 . One may verify (I3.3p . remarkably using the 
evolution (|3.10p alone. We also have the gauge factor, given by 

u n = - hb 2 t 2 (w° n+1 - w n ) 2 (w° - w n+1 ) -f- bit (w° +1 - u£) 

(a 2 (w° +1 - wl l+1 ) (w n+1 - wl) + b 2 t (w° +1 - wi) (w° - w n+1 )) 
- a i02 (w n - w n+1 ) (w n - w* +1 ) 2 J , 
which satisfies the equation 

w„ +2 [z n - l){0i{z n {a 2 y n - 1) + 1) - a 2 tbib 2 y n ) 



i) 

u n y n (0iz„ - tbib 2 )(a\y n z n + 9ity n - a 2 (t(bi + b 2 )y n + z n - 1)) ' 
We may now parameterize the deformation matrix for the double shift in terms of 



y n , z n and w n as 



s(fli-tfeife 2 ) 



r>( +\-TJ-\ b 1 +b 2 + (C n -y n )e 1 b 1 +b 2 + ((; n -y n )8 1 

tS(X, t) - I + I x((tb 1 fa tl -c„)-i)9i-^)((t&2(y-^)-i)9i-«^) ^162(61 +h2+(Cn-y re )fli)t 2 



<5„w„(ei-t6i6 2 )82 6\—tbxb 2 

Curiously, the coefficient of x is lower triangular and the constant coefficient is 
upper-triangular. This rather simplified Lax pair comes at the expense of the 
requirement that the variables y n and z n explicitly lie on the biquadratic 

V(y, z, t) =9i(z - i f - 9 iy (z - 1) ((oi + a 2 ) z - t {h + b 2 )) 

+ y 2 (aia 2 z - 6 it) (9iz - tbib 2 ) = 0, 

which makes the computation of the compatibility condition, 

A(x,q 2 t)B(x,t) = B(q 2 x,t)A(x,t), 

slightly more difficult. One can verify directly that V(y,z,t) = implies that 
V(y', q 2 z, q 2 t) = and V(y', z', q 2 t) = under the evolution defined by ([3~3]) , This 
forms an explicit parameterization of the member of the pencil of biquadratics for 
each t. The existence of such a parameterization is not without precedent, and has 
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appeared in the work of Yamada [41] and Noumi et al. [42] . We do remark that 
it is interesting that this explicit dependence on the bi-quadratic curve essentially 
came from a condition on the Lax matrices. 

3.2. q-P(A^) as a reduction of dSKdV. While a reduction to q-P(A { 3 1] ) 
has been provided by one of the authors previously [30] , we know of no reduction 
from an integrable lattice equation to any member of the hierarchy above q-P(A^). 
We wish to extend this further and consider reductions of a Painleve equation that 
is higher up in the classification scheme; namely the q-Painleve equation associ- 
ated with the A 2 1} surface, (or q-P(A^)), given by We will find that this 
equation appears as a reduction of (|3.1[) , or Ql^°, or the non-autonomous discrete 
Schwarzian Korteweg-de Vries equation [43] . 
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Figure 5. The reduction and the labelling of variables. 

We impose (]0.2[) on (I3.1[) with s\ = 4 and s 2 — 2, hence, our constaint, (jl.llbl) . 
becomes 

O-l+i Pm+2 4 

= —3 — : = 1 . 

Oil p m 

which we solve in a similar manner as before; 



a\q l if I — mod 4 

a 2 q l if 1 = 1 mod 4 __ f b x q 2m if m = mod 2 

a 3 q l if I = 2 mod 4 ' 1 m ~ \ b 2 q 2m if m = 1 mod 2 

a^q 1 if I = 3 mod 4 



(3.13) ai 



We use the same Lax pair for (]3.1|) as the previous section, namely (|1.4|) where ii, m 
and Mj jm are specified by (|3.2|) . The i-direction is also chosen to be constant in 
the direction of the reduction. We choose our spectral parameter, x, and indepen- 
dent variable, t, in a manner in which correspondence with g-P^j 1 '') comes more 
naturally, that is 

t = q 2m - l =q n , x = ( t = t. 

1 1 

We note that the generating shift is no longer just a shift in m alone, but a si- 
multaneous shift in I and m , i.e., (I, m) (J, + 1, m + 1), hence we present a Lax 
pair that represents the shift (l,m) — > (I + 4, m + 2), which is now equivalent to 
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(x,t) — > (q 4 x,t), and the shift (l,m) — > (I + l,m+ 1), which is equivalent to the 
shift (x,t) —> {qx,qt). Hence, our Lax pair is a linear system satisfying 

(3.14a) Y(q 4 x, t) = A(x, t)Y{x, t), 

(3.14b) Y(qx, qt) = B(x, t)Y(x, t), 

where we compose A and B in terms of L and M in the following way 

(3.15a) 

A(x, t) ■ ( r-lLi +3 m+ 2Li + 2^ m +2Ml+2,m+lLl + l^n+lLL m +lM^ m , 

(3.15b) A{x,t) = 




u n+2 w n+l 

The deformation matrix may be written as 

(3.15c) B{x,t) ^ L Um+1 M Lm , 

(3.15(1) 



1 W° n +2-<+l\ I 1 <-<+2* 



u n+2 w n+l 

B(x,t) = I xa i i 

1 n i -L 



n+2 "'ra+l 




1 



where we have written these as functions of x and t. The compatibility of the 
system defined by (|3.14[) is given by 

(3.16) A (qx, qt) B(x, t) = B(q 4 x, t)A(x, t). 

One may use (|3.16l) to derive the required evolution equations; the equations defin- 
ing the evolution may be written as 

/o 17 ^ _ q 4 «-2 ~ W n) + b 2 tw n ( W n-1 " K-2) 

(o.iia) u ' n +i — 1 n a , , ,/ 1 o^ ' 

«2<-l (w„-2 ~ W°) + ht - W^_ 2 ) W° 

w x n _2 (a 2 - bit) - a 2 w° + Mw°_i 

. . (2 2 Q3 CZ4 o 

(3.17c) ai — > — , a 2 — > — , 03 — > — , 04 — > a\q , 

9 9 9 

(3.17d) 6i->-4 ) b 2 ^q 2 b 1 . 

q 2 

The difficult step is to extract a second order system from this seemingly sixth order 
system via some special parameterization. To do this, we observe the properties 
of (|3.14l) as an associated linear problem. Firstly, expanding out (13.15b[) in the 
spectral parameter, we find 

A(x, t) = A + A x x + A 2 x 2 + A 3 x 3 , 



(3.17b) w 



n+l 
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where Ao is the identity matrix and A3 is a lower triangular matrix with diagonal 
entries 9\t and # 2 i, where 

(< - W r\+l) (<+2 ~ O ( W n+1 ~ W n+ 2 ) 

_ aia 4 6 2 « - w„+i) « - K+2) K+i - ^+2) 



where, 6\ and 62 are invariants in accordance with the evolution equations, given 
by (|3.17p . The simplicity of the individual factors of A(x,t) from (|3.15bl) give us 
that the determinant is 

(3.18) det A(x, t) = (xa t - 1) (xa 2 ~ 1) (xa 3 - 1) (xa 4 - 1) (txbi - 1) (txb 2 - 1) . 

These properties should remind us of the properties of a special case of the spectral 
matrix in the Lax pair of Sakai 1 3 5 j . We may simply use a transformation of the 
form Y(x, t) — > SY(x, t), where S is a lower triangular matrix that is constant in x 
such that the transformation A — > SAS^ 1 diagonalises A3. While this matrix, S, 
isn't particularly nice to write down, the resulting matrix is in the general form 

(Bit {{x - y n )(x - e n ) + C,n) 9 2 tuj n (x-y n ) \ 
O^jx-Jn + 6 n ) 
02t{(X -y n ){x - Xn) +Vn) 
w n / 

The terms e n , Xm In and 5 n may be determined from (|3.18j) in terms of y n and Cm 
r] n , the Oi's and the b^s. The relation between this form of A(x,t) and the known 
Lax pairs of Sakai [35] and Yamada |41j has recently been found by one of the 
authors |44j . Following the motivation from the previous section, this parameter- 
ization defines y n to be the root of the (l,2)-entry, leaving a choice of £ n and rj n 
such that 

{OlKnVn + l)(Q2tT) n y n + 1) = dct MVn , t) , 

where the right hand side is defined by (|3.18|) . With this in mind the defining 
equations for C, n and r\ n are 

OitCnVn + 1 = (1 - y n z n ) (hty n - 1) (b 2 ty n - 1) , 

a . . , (aiVn - 1) (a 2 y n - 1) {a 3 y n - 1) («4j/n - 1) 

V2tr] n y n + 1 = r . 

(1 - y n z n ) 

This specifies y n and z n in terms of the lattice variables: 
(3.19) 



Vn 



w l+2 - wi +1 )((w° - - w r \ +2 )(fe 2 tK° - - «^+ 2 ) 

- a 3 (w° +1 - wl)(w® - u£ +2 )) + a 2 {w° n - w*)(wl - w* +1 )(w° - u£ +2 ) 

( W «+2 - W n+l)) + a l( W n ~ W «+l)( W n ~ W li+l)( W n+2 ~ W n)( W n ~ W n+2) 

(w° n+1 - w l n+2 ) ai(w° - wl l+1 )(w° +2 - u£)(a 3 (u£ +1 - u£) 

K° - w I 1 1+2 )(w^ +2 - u£+i) + - wl+i)(wl - w n+2 ) 

(w° l+1 - wl l+2 )) - a 2 a 3 (w° +1 - wl)(w* - u>i + i)(«£ +2 - 

K-^ +2 ) 21 
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(3.20) 

z n =[ai(w° n - w^ +1 )(K° l+2 - wl +2 )(a 3 {w° +1 - w x ^){uP n - i4 +2 )(w* +2 - u£+i) 
+ ht(w n - - w, 1 l+2 )(w° +1 - w] l+2 )) - a 2 {w° n+2 - w l n+1 ) 

(w ( n - wl l+2 )(wl l+2 - w l n ){w l n+2 - + a 2 a 3 (w° +1 - w*)(w£ +2 _ W n+l) 

{w l n+1 - w] l+2 ){wl - w^ +2 ) 2 ] [{w° n - w l n+2 ){w l n+2 - w n+l ){{w° n - w] l+1 ) 
{ai{w° n+2 - wDiw 1 ^ - w\ +2 ) - b 2 t(wl +2 - - wl +2 )) 

+ a 2 (wl l - w,\ + i)(w° +2 - wl +1 ){w° - wl +2 ))], 

which satisfy (|0.3a|l and (|0.3b[) under the identification of y = y n , z = z n , ?/ = Vn+i 
and z' = z„+4. While these expressions may fail to be succinct, they do succeed in 
being very explicit. Furthermore, using this identification, it is possible to directly 
verify (Iu3a1) and ((uT3b| from ([3~P7| alone. 

4. Conclusion 

Given a partial difference equation of the ABS-list, we have outlined a way 
in which we may find a Lax representation for any similarity (travelling wave) 
reduction of the form (j0.2|) . The method outlines how the Lax matrices may be 
expressed in terms of products of the Lax matrices of the partial difference equation 
from which it was derived. 

We have broken new ground by deautonomizing this theory so that we were 
able to give an algorithmic way to derive Lax representations of reductions of deau- 
tonomized lattice equations. We have also found the relation between the derived 
reductions and known Painleve equations. 
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